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Abstract: A new non-Hermitian E2-quasi-exactly solvable model is constructed con¬ 
taining two previously known models of this type as limits in one of its three parameters. 
We identify the optimal finite approximation to the double scaling limit to the complex 
Mathieu Hamiltonian. A detailed analysis of the vicinity of the exceptional points in the 
parameter space is provided by discussing the branch cut structures responsible for the 
chirality when exceptional points are surrounded and the structure of the corresponding 
energy eigenvalue loops stretching over several Riemann sheets. We compute the Stielt- 
jes measure and momentum functionals for the coefficient functions that are univariate 
weakly orthogonal polynomials in the energy obeying three-term recurrence relations. 


1. Introduction 


In addition to the interesting mathematical aspect of enlarging the set of sl 2 {C) 01 to 
£^ 2 -quasi-exactly solvable models Q, the latter type also constitutes the natural framework 
for various physical applications in optics where the formal analogy between the Helmholtz 
equation and the Schrodinger equation is exploited i I, i 0 g 10,0, in, 0 • Fur¬ 
thermore, a special case of these systems with a specific representation corresponds to the 
complex Mathieu equation that finds an interesting application in nonequilibrium statisti¬ 
cal mechanics, where it corresponds to the eigenvalue equation for the collision operator in 


a two-dimensional classical Lorentz gas [14 


Here we are mainly concerned with the extension of quasi-exactly solvable models 
11,0,11, 0 , g to non-Hermitian quantum mechanical systems |^, 0,0, [0 within the 
above mentioned scheme. So far two different types of S 2 -uiodels have been constructed 


in jg, 241 and the main purpose of this manuscript is to investigate whether it is possible 
to construct a more general model that unifies the two. We show that this is indeed 
possible by combining the two models and introducing a new parameter into the system 
that interpolates between the two. In a similar fashion as the previously constructed 
models, also this one reduces in the double scaling limit to the complex Mathieu equation. 
As that equation is not fully explored analytically this limit provides an important option 
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to obtain interesting information about the complex Mathieu system. On the other hand, 
for some applications it may also be sufficient to study an approximate behaviour for some 
finite values of the coupling constants. For that purpose we identify the parameter for 
which the general model is the optimal approximation for the complex Mathieu system. 

Our manuscript is organized as follows: In section 2 we introduce the general unifying 
model involving three parameters. We determine the eigenfunctions by solving the stan¬ 
dard three-term recurrence relations for the coefficient functions and determine the energy 
eigenfunction from the requirement that the three-term recurrence relations reduce to a 
two-term relation. We devote section three to the study of the exceptional points and their 
vicinities in the parameter space. The explicit branch cut structure is provided that ex¬ 
plains the so-called energy eigenvalue loops. In section 4 we compute the central properties 
of the weakly orthogonal polynomials entering as coefficient functions in the Ansatz for the 
eigenfunctions, i.e. their norms, the corresponding Stieltjes measure and the momentum 
functionals. We state our conclusions in section 5. 


2. A unifying E2-quasi-exactly solvable model 


The general notion ii underlying solvable Hamiltonian systems is that its Hamiltonian 
operators % acting on some graded space Vn as, % : Vn Vn preserves the flag structure 
Vo C Vi C F 2 C ... C Fn C ... A distinction is usually made between exactly and quasi- 
exactly solvable, depending on whether the structure preservation holds for an inhnite or 
a finite flag, respectively. Here we are concerned with the latter. Lie algebraic versions of 
Hamiltonians in this context are usually taken to be of s/ 2 (C)-type [||, but as recently 
proposed |P, they may also be taken to be of a Euclidean Lie algebraic type, thus 
giving rise to qualitatively new structures. 

At present two different types of E 2 -quasi-exactly solvable models were identified 


Hg = j2 + C2(u2-u2)2 + 2iCA(u2-u2), (2.1) 

+ CuvJ + 2iCN{u^ - v^), (2.2) 


in Q and |^|, respectively. Both Hamiltonians are expressed in terms of the E 2 -basis 
operators u, v and J that obey the commutation relations 


[u, J] = iv, [v, J] = —iu, [u, v\ = 0. 


(2.3) 


Except for at N = 1/4, both Hamiltonians are non-Hermitian, but respect the anti- 
linear symmetry |2^ VT 3 : J ^ J, u ^ v, v ^ u, i ^ —i as defined in |^|. Eor the 
particular representation J := —idg, u := sin0 v := cos 6 the PTa-symmetry is simply 
VT 3 : 0 —>■ 7r/2 — 9, i ^ —i, such that the invariant vector spaces over R were dehned as 


Vn{4>o) ■ = span { (pQ [sin(20), isin(40),..., sin(2n0)] | 9 G R,'PT 3 (</>o) = (/)o € L} , (2.4) 

^n{4>o) '■ = span{(/)Q [l,icos(20),...,i”cos(2n0)]| 9 € R,PT 3 (())o) = (pQ ^ L} ■ (2-5) 


In order to construct Hamiltonians that preserve the flag structure one needs to identify 
the action of the E' 2 -basis operators and its combinations on these spaces as explained in 
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more detail in |Q. The behaviour found allowed to identify the Hamiltonians and 


^£2 ™ (P-lD and as quasi-exactly solvable. The general structure suggests that there 
might be a master Hamiltonian that unifies the above Hamiltonians into one preserving the 
quasi-exact solvability. We demonstrate here that this is possible and study the properties 
of that model. 

Thus we introduce the new Hamiltonian 


C, A) = + 2(1 - X)CuvJ + XC {u^ - v^Y + 2 iCN{u^ - X,C,N e M, (2.6) 


and demonstrate explicitly that it is indeed ill 2 -quasi-exactly solvable. First we observe 
that T-LiN, A) interpolates between the two models in (^) and (|2.2|) by varying A, since 


hmn{N,c,x) = n^^l 


and 


lim'H(2iV,C/2,A)=il^g. 

A—^0 ^ 


(2.7) 


Furthermore, 77(iV, A) reduces to the complex Mathieu Hamiltonian in the double scaling 
limit lim7v-).oo,c^o 77(A^, C, A) = T^Mat = + 2 ig{u^ — u^) for g := NC, < 00 . We also note 

that 77^(A^, C, A) = 77(1 —A —iV, A), which implies that 77(iV, C, A) is non-Hermitian unless 
2N = 1 — A, with free coupling constant G M. 

Given the structure for the vector spaces in (^^) and (^]^) we now make the follow¬ 
ing Ansatze for the two fundamental solutions of the corresponding Schrodinger equation 


= ^ 0 '^ i'^CnPn{E) cos{2n0), and = ^ 0 '^ i^'^^CnQniE) sm{2n0), (2.8) 

n=0 n=0 


where the T’Ts-symmetric ground state is taken to be (/>q = and the constant 

Cn is Cn = 1/C^{N + A)(l -|- A)”' ^ [(1 -t- + 2A)/(1 + A)]^_]^ with (a)n := F (a -|- n) /F (a) 

denoting the Pochhammer symbol. The constants are chosen conveniently in order to 
ensure the simplicity of the to be determined n-th and (n — l)-th order polynomials Pn{E), 
Qn{E) in the energies E, respectively. Upon substitution into the Schrodinger equation we 
obtain the three-term recurrence relations 


P2 = {E- X(Y - 4)Pi + 2Y [N - 1] [A^ + A] Po, (2.9) 

Pn+i = {E- \Y - 4n2)P„ + [A^ + nA + (n - 1)] [N - {n - 1)A - n] P„_i, (2.10) 

Q2 = {E-^-Xe)Qi, ( 2 . 11 ) 

Qm+I = {E - XY - 4m^)Qm + [A^ + mX {m - 1)] [A^ - (m - 1)A - m] Qm-i,(2.12) 

for n = 0, 2,... and for m = 2,3,4,... Note that a more generic Ansatz for the unifying 

model involving two independent coupling constants /i, A in the terms iJiQuvJ-\-X(Y{u^—v‘^Y 
leads to a four term recurrence relation in which the highest term is always proportional to 
// -|- 2A — 2. Thus taking this term to zero with the appropriate choice for // reduces this to 
the desired three term relations that may be solved in complete generality as outlined in 
Q. The lowest order polynomials are easily computed in a recursive way. Taking Pq = 1 
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we obtain 


Pi= E- AC^ (2.13) 

P2 = + 2C^ [A - AP + N{\ + N -1)] + {E- A)E, 

P3 = -A^C^ + (a(2A + 3P - 13) - - 3(A - l)iV + 2) + (P - 16) (P - 4)P 

-C^ [3Ap2 + P (2A2 - SiV^ - 3A(Ar + n) + 3iv + 2) + 32(A + N{\ + N - 1))] , 
and likewise with Qi = 1 we compute 

Q2 = E-A-XC\ (2.14) 

Qz = [A(15 - 2A - 2P) + N‘^ + {\- 1)A^ - 2] + (P - 16)(P - 4), 

= -X^cf^ + AC^ [8 + A(8A + 3P - 38) - 2N‘^ - 2(A - l)iV] + (P - 36)(P - 16)(P - 4) 

+C^ [-8 {-12X^ + 69A + 5AiV + 5(iV - l)N - 12)] 

+C^ [-3AP2 + 2P ((47 - 4A)A + + {X - 1)N - A)] . 


In both cases we observe the typical feature for quasi-exactly solvable systems that the 
three term relation can be reset to a two-term relation at a certain level. This is due to the 
fact that in ( 2.1C| ) and ( |2.12| ) the last term vanishes when m = n = h = —(l + iV)/(l + A) 
or m = n = h= (A + N) /(I + A). Thus when taking N = n-\- (n — 1)A we find the typical 
factorization 


and (2.15) 

The first solutions for the factor are easily found from ( |2.1C| ) and ( |2.12| ) to 

Pl = P-4h2-AC^ (2.16) 

R2 = {E- 4n^ - AC^)(P - 4(h + 1)^ - AC^) - 2h(l + A)^^^. (2.17) 


Next we compute the energy eigenvalues En from the constraints Ph{E) = 0 and 
Qn{E) = 0 for the lowest values of N. For the solutions related to the even fundamental 
solution in (^) we find 

N = 1: El = AC^ 


iV = 2 + A ; 
iV = 3 + 2A 


Po’=^ = 2 


+ AC' ± 2^Jl-{l + XYf\ 


Ef = f + AC^ + f + 1 [52 - 12(1 + A)^Cl A'S 


(2.18) 

(2.19) 

( 2 . 20 ) 


with 1)3 := 35 -b 18(A + [3(A -b - 13] ^ -b [18(A + + 35] ^ £ = 0, ±2. 

For the solutions related to the odd fundamental solution in ( p.8D we obtain 

P = 2-bA: P| = 4 + AC', (2.21) 

iV = 3 + 2A: P^’± = 10 + C'A ± 2^9 - (A + l^C^, (2.22) 

iV = 4 + 3A: Ef = — + XC'^ + —e^+ -[l96-l2{l + XfC^]e-^n-^, (2.23) 

3 3 3 

with fl3 := 143 -b 18C'(A -b 1)^ + (3C'(A + 1)2 - 49)^ + (l8C'(A + 1)2 + 143)^ i = 0, ±2. 

Solutions for higher order may of course also be obtained, but are rather lengthy and 
therefore not reported here. 
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3. Exceptional points and their vicinities 


The special point in parameter space where two real energy eigenvalues viewed as functions 
of the coupling constants merge and subsequently split into a complex conjugate pair is 
usually referred to as exceptional point |2^, ^|. In our system these points can be 

computed in an explicit simple and straightforward manner. Using that by definition the 
discriminant A equals the product of the squares of the differences of all energy eigenvalues 
for 1 < i < n, i.e. A = ni<i<j<n(-^* “ obtains the exceptional points from 

the real zeros of A(£'). For practical purposes one may also exploit the fact P], that 
the discriminant equals the determinant of the Sylvester matrix. This viewpoint has the 
advantage that it does not require the computation of all the eigenvalues and is more 
efficient when the sole purpose is to find the exceptional points. Thus in our case we have 
to find the real zeros of the discriminants A? and A? for the polynomials Ph{E) and Qfi{E), 
respectively. Extracting overall constant factors k as A = kA, that do not contribute to 
the zeros, we obtain for the lowest values of h 


Ai = ? - 1, (3.1) 

AI = c'-9, 

A^ = C" - C" + lOsf - 36, 

A^ = - 37^"^ + 991C^ - 3600, 

+ 385C® - 33120C^ + 16128C^ - 732276C^ + 129600, 

A^ = - 94C^° + 7041C® - 381600^*^ + 6645600^"^ - 78318900C^ + 158760000, 

where we abbreviated C := C(1 + ^)' 

There exist many detailed studies about the structures in the coupling constant space 

, Q . It is evident that when tracing 

a complex energy eigenvalue E as functions of the coupling constants, A or in our case, 
the corresponding path in the energy plane will inevitably pass through various Riemann 
sheets due to the branch cut structure. As a consequence one naturally generates eigenvalue 
loops that stretch over several Riemann sheets. This phenomenon is well studied for a large 
number of models and we demonstrate here that it also occurs in quasi-exactly solvable 
models. The basic principle can be demonstrated with the square root singularity occurring 
in with branch cuts from (— oo,— 1 — 1/C) and (1/C — l,oo). The energy loops are 
generated by computing = A + pe^'^‘^,C) for some fixed values of C, center A and 

the radius p in the A-plane as functions of 4> as illustrated in hgure |^a) and (b). In panel 
(a) we simply trace the energy around a point in parameter space that leads to two real 
eigenvalues. For a small radius ones reaches the starting point by encircling A just once. 
However, when the radius is increased one needs to surround A twice to reach the starting 
point and when the radius is increased even further one only needs to surround A once 
switching, however, between both energy eigenvalues. 

Essentially this structure survives when the two eigenvalues merge into an exceptional 
point. However, since the exceptional point is a branch point we no longer have the option 
for a closed loop around it produced from only one energy eigenvalue as seen in figure ||(b) . 


in the vicinity of the exceptional points 


31 
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Re(E) 


Re(E) 


Figure 1: Energy eigenvalue loops E 2 ^{\ + around two real eigenvalues panel (a) and 

around an exceptional point panel (b) as functions of indicated by the numbers on the loops, for 
fixed value of C = 1/2 at A = 1/10 in (a) and A = 1 in (b). The energy eigenvalues for p = 0 are 
distinct in panel (a) as = 0.35, = 3.70 and coalesce to an exceptional point in panel (b) 

as ^2 ■ = ^ 2 ’+ = 9/4. 
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Figure 2: Energy levels and branch cut structure for E^’^ for fixed (/ = 1/2 as functions of A. The 
branch cuts extend to the left and right from the exceptional points (—oo, —3) and (1, oo). 


This behaviour is easily understood from the structure of the branch cuts as depicted 
in figure]^. Whereas for small radii it is possible to encircle for instance the point A = 1/10 
without crossing any branch cut, this is not possible when encircling the exceptional point 
at A = 1 where we have to analytically continue from to ^ 2 ^ when crossing a cut. 
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This structure is the same for intermediate radii. For large radii we cross the first cut 
already at a half circle turn, such that one returns back to the original value already after 
one complete turn. 

When more eigenvalues are present the structure will be more intricate. Considering for 
instance a scenario with four eigenvalues in the form of two complex conjugate eigenvalues 
and an exceptional point, see figure |3|(a), we need to perform again at least two turns in 
the A-plane in order to return to the initial position for the energy loops when surrounding 
an exceptional point. The two complex conjugate eigenvalues may be enclosed with just 
one turn, albeit we require again different energy eigenvalues for this. When enlarging 
the radius the loops will eventually merge as depicted in figure |^(b) for a situation with 
a degenerate complex eigenvalue and two complex eigenvalues. We observe that for the 
given values we have to surround the chosen point at least three times to obtain a closed 
energy loop surrounding the indicated centers. 




Re(E) Re(E) 

Figure 3: Energy eigenvalues El{X + Q) as functions of </>, indicated by the numbers on the 

loops, for fixed value ( = 1/2 at A = 9.5284 in (a) and A = 5.2562 + 19.9526 in (b). The energy 
eigenvalues for p = 0 in panel (a) are = Ef^ = 25.6613, Ef^ = (Ef‘^)* = 7.1029 + 129.8106 
and Ef^ = Ef^ = 37.7449-18.7611, Ef^ = 9.8103 + 16.7668, Ef"^ = -24.0439 + 120.7081 in panel 
(b). The radii are p = 4.0 and p = 8.5 in panel (a) and (b), respectively. 


In the same manner as for the simpler scenario one may understand the nature of these 
loops from an analysis of the branch cut structure of the energy as seen in figure Tracing 
the indicated radii at p = 4.0 and p = 8.5 in figure ^ produces the energy loops in hgure ^ 
when properly taking care of the analytic continuation at the branch cuts. 


As discussed earlier the Hamiltonian T-L{N, (, A) has the interesting property that in the 
double scaling limit it reduces to the complex Mathieu equation for which only incomplete 
information is available, especially concerning the locations of the exceptional points. In 
comparison with the previously analyzed models in Q and in |24] we have now the 
additional parameter A at our disposal and we may investigate how the complex Mathieu 
system is approached. In particular we may address the question of whether there exists a 
value A for which this is optimal. Our numerical results are depicted in figure We find 
a similar qualitative behaviour for the other exceptional points, which we do not report 
here. 
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Comparing the rate of the approach for different values of A we conclude that A 

1) is the best approximation to the complex Mathieu system for some finite values of N. 
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If one is exclusively interested in the computation of the exceptional point it is most 


efficient to carry out the double scaling limit already for the three-term relation (2.10) and 
( 2.12 ) as explained in ^|. 



n 


Figure 5: Double scaling limit of limAr->.oo,c ->-0 A) = 'Huat to the smallest exceptional point 

at = 1.46877 with A(n) = Co-^(^) ~ Cmj N{n) = (n -I- 1) -I- nA for n = 1, 2, 3,_ 


4. Weakly orthogonal polynomials 

It is well known from Favard’s theorem ||3^ , ^ that polynomials ^n{E) constructed from 
three-term relations in the way mentioned above possess a norm 

= N^dnm. (4.1) 


defined by the action of a linear functional L acting on arbitrary polynomials p in as 

poo 

Cip) = / p{E)uiE)dE, C{1) = 1. (4.2) 


This norm may be computed in two alternative ways. The simplest way is to multiply 
the three-term relation by 4>„_i and act subsequently on the resulting equation with C. 
Using the property together with (O) then simply yields 

UL -i&fc) where the are the negative coefficients in front of Whereas the 

first method simply assumes that the functional exist the second method goes further and 
actually provides an explicit expressions for the measure. As argued in the concrete 
formulae for uj{E) may be computed from 


= '^u}k6{E - Ek 


(4.3) 


k=l 
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where the energies are the (. roots of the polynomial ^{E). The ^ constants ojk can be 
determined by the i equations 

i 

'^uJk^n{Ek) = 5no, forn E No- (4.4) 

fc=i 

In our case the integer are determined from N = — 1)\ and N = {(. -\- 1) -\- (.X for 

the Pi{E) and Qi^i{E), respectively. 

Using the first method we obtain 


= 2C^^{1 + Xf^ 
1 


n9 = 


l-N 
1 + A 


-Nf 


X + N 
1 + A 


n = 1, 2, 3,... 
n = 2,3, 4, 


(4.5) 

(4.6) 


2(iV +A)(l - A^)' 

with = A^^ = 1. Due to the non-Hermitian nature of the Hamiltonian this norm is in 
general not positive definite. For instance for N = 4 + 3A we have 

N[ = 1, N[ =-24C^{1 +Xf, A^f = 240C^1 + A)^, A^l" =-1440C®(1 + A)®. (4.7) 

The exception is the class of models where the Hamiltonian becomes Hermitian, i.e. when 
A = 1 — 2N holds. For this value of A the expressions in (|4.5|) and (^) become positive 
definite 


nP ^ = 2C^(A^ - ifN^. 


(4.8) 


Let us now consider the second method and compute explicitly the measure for a few 
examples. For A^ = 2 + A and N = 3 + 2A we solve (^) for the even and odd solutions, 
respectively, to 


a;± = - ± — =, and = - ± — ■ 

2^1-(1 + A)2c2 2^9-(1 + A)2c2 

Computing now (|4.1| ) with (|4.2D agrees with (4.5) and (|4.6D 

nP = C{P^) =uj\ + ujf = l 


N( = C{Pl) = E^f^ - AC" + cui E^f- - AC" = -4C , 


= C{Qi)=uj%[El^^-4-XC] +ujf[El^--4-XC] =-AQ 


(4.9) 

(4.10) 

(4.11) 

(4.12) 


Similarly we compute for N = 3 + 2A 




260 - 60C^) H + (^3C^ + 4 ) H2 + 20H3 


12 


3 + 


13 - 3^^)^ + fl3 - 3? ) H2 + 


+ 


0 J 2 = X-2> ^3 = X2> (4-13) 


^ (3C -20H + 4] (l + 2e- 


■7.2 iTri ^ 

4 + 3C -20e—H 


36(3C + H2 - 13) 


12 1 + 2e~ 


-13 + 1-e^ H2 
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and confirm that 


Nq — C{Pq) — + ^2 + W3 — 1, 

< = C{Pl) = ulPl{E'f^)+ojlPKEr^)+ojlPKEf) = -12C\ 


(4.14) 


< = E{Pi) = ulPhE'f^) + ojlP^Er^) + ojlPhEf) = 48C 
C{P^P2) = u:\P^{Ef)P^{Ef)+ulP^{El^-^)P^{E'f-^)+ulP^{^^^^^ 


= 0 . 


Note that the last relation in (4.14) does not follow from the first method. 


As the final quantity we also compute the moment functionals defined in |^, as 


1^1 


-.= C(E") 


n—1 


(”) , 


k=l 


= l^k, 

k=0 


(4.15) 


Once again also these quantities can be obtained in two alternative ways, that is either 
from the computation of the integrals or directly from the original polynomials Pn and Qn 

(n) 

without the knowledge of the constants ca^. In the last equation the coefficients ^ are 
defined through the expansion Pn{E) = 2'^~^E^ — X)fc=o E^ and Qn{E) = — 

X)fc=o for our even and odd solutions, respectively. For the even solutions with 

N = 2 + A we obtain 


1^0 = 1, (4.16) 

/if = ACf (4.17) 

/if = AV-4f, (4.18) 

/if = A^C® - 12AC2 c' - lef, (4.19) 

/if = A^C® - 24X^C^C‘^ + 16 (C^ - 1) ^ - 64^, (4.20) 

and similarly for the odd solutions with N = 3 + 2A we compute for instance 

/i® = 1, (4.21) 

^f = 4 + ACf (4.22) 

/if = 16 - 4? + A^Cf (4.23) 

/if = A^C® - 12 (A^ + A^ + A) - 48(2A2 + 3A + 2)C^ + 64. (4.24) 


Thus P{N, (, A) possesses indeed all the standard features of a quasi-exactly solvable 
model of i72-type. 

5. Conclusions 

Following the principles outlined in |^] we have constructed a new three-parameter quasi- 
exactly solvable model of i72-type. One of the parameters can be employed to interpolate 
between two previously constructed models. With regard to one of the original motivations 
that triggered the investigation of these models, that is the double scaling limit towards 
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the complex Mathieu equation, we found that for A = 1, i.e. finite values for N 

best approximate the complex Mathieu system and mimic its qualitative behaviour. We 
provided a detailed discussion of the determination of the exceptional points and the en¬ 
ergy branch cut structure responsible for the intricate energy loop structure stretching 
over several Riemann sheets. The coefficient functions are shown to possess the standard 
properties of weakly orthogonal polynomials. 

Acknowledgments: I am grateful to Kazuki Kanki for making reference |T^ available to 
me. 

References 

[1] A. V. Turbiner, Quasi-Exactly-Solvable problems and sl(2) Algebra, Commun. Math. Phys. 
118, 467-474 (1988). 

[2] A. Turbiner, Lie algebras and linear operators with invariant subspaces, Lie Algebras, 
Cohomologies and New Findings in Quantum Mechanics, Contemp. Math. AMS, (eds N. 
Kamran and P.J. Olver) 160, 263-310 (1994). 

[3] A. Fring, F2-quasi-exact solvability for non-Hermitian models, J. Phys. A48, 145301(19) 
(2015). 

[4] Z. H. Musslimani, K. G. Makris, R. El-Ganainy, and D. N. Christodoulides, Optical Solitons 
in PT Periodic Potentials, Phys. Rev. Lett. 100, 030402 (2008). 

[5] K. G. Makris, R. El-Ganainy, D. N. Ghristodoulides, and Z. H. Musslimani, PT-symmetric 
optical lattices, Phys. Rev. A81, 063807(10) (2010). 

[6] A. Guo, G. J. Salamo, D. Duchesne, R. Morandotti, M. Volatier-Ravat, V. Aimez, G. A. 
Siviloglou, and D. Ghristodoulides, Observation of PT-Symmetry Breaking in Complex 
Optical Potentials, Phys. Rev. Lett. 103, 093902(4) (2009). 

[7] B. Midya, B. Roy, and R. Roychoudhury, A note on the PT invariant potential 
Acos'^x ‘iiVQSin2x, Phys. Lett. A374, 2605-2607 (2010). 

[8] H. Jones, Use of equivalent Hermitian Hamiltonian for PT-symmetric sinusoidal optical 
lattices, J. Phys. A44, 345302 (2011). 

[9] E. Graefe and H. Jones, PT-symmetric sinusoidal optical lattices at the symmetry-breaking 
threshold, Phys. Rev. A84, 013818(8) (2011). 

[10] S. Dey, A. Fring, and T. Mathanaranjan, Non-Hermitian systems of Euclidean Lie algebraic 
type with real eigenvalue spectra. Annals of Physics 346, 28-41 (2014). 

[11] S. Dey, A. Fring, and T. Mathanaranjan, Spontaneous PT-symmetry breaking for systems of 
noncommutative Euclidean Lie algebraic type, arXiv: 1407.8097, to appear in Int. J. Theor. 
Phys. 

[12] S. Longhi and G. Della Valle, Invisible defects in complex crystals. Annals of Physics 334, 
35-46 (2013). 

[13] K. Makris, C. D. Musslimani, Z.H., and S. Rotter, Constant-intensity waves and their 
modulation instability in nonHermitian potentials, arXiv:1503.08986. 


- 12 - 



A unifying E2-quasi-exactly solvable model 


[14] K. Kanki, Spontaneous breaking of a PT-symmetry in the Liouvillian dynamics at a 
nonhermitian degeneracy point, talk at the 15th International Workshop on 
Pseudo-Hermitian Hamiltonians in Quantum Physics, May 18-23, University of Palermo, 

Italy (2015). 

[15] Z. Zhang, Irreversibility and extended formulation of classical and quantum nonintegrable 
dynamics, PhD Thesis, The University of Texas at Austin (1995). 

[16] A. Khare and B. P. Mandal, A PT-invariant potential with complex QES eigenvalues, Phys. 
Lett. A 272, 53-56 (2000). 

[17] B. Bagchi, S. Mallik, C. Quesne, and R. Roychoudhury, A PT-symmetric QES partner to the 
Khare-Mandal potential with real eigenvalues, Phys. Lett. A 289, 34-38 (2001). 

[18] C. M. Bender and M. Monou, New quasi-exactly solvable sextic polynomial potentials, J. 
Phys. A 38, 2179-2187 (2005). 

[19] B. Bagchi, C. Quesne, and R. Roychoudhury, A complex periodic QES potential and 
exceptional points, J. Phys. A 41, 022001 (2008). 

[20] F. G. Scholtz, H. B. Geyer, and F. Hahne, Quasi-Hermitian Operators in Quantum 
Mechanics and the Variational Principle, Ann. Phys. 213, 74-101 (1992). 

[21] C. M. Bender and S. Boettcher, Real Spectra in Non-Hermitian Hamiltonians Having PT 
Symmetry, Phys. Rev. Lett. 80, 5243-5246 (1998). 

[22] C. M. Bender, Making sense of non-Hermitian Hamiltonians, Rept. Prog. Phys. 70, 947-1018 
(2007). 

[23] A. Mostafazadeh, Pseudo-Hermitian Representation of Quantum Mechanics, Int. J. Geom. 
Meth. Mod. Phys. 7, 1191-1306 (2010). 

[24] A. Fring, A new non-Hermitian E2-quasi-exactly solvable model, Phys. Lett. 379, 873-876 
(2015). 

[25] E. Wigner, Normal form of antiunitary operators, J. Math. Phys. 1, 409-413 (1960). 

[26] T. Kato, Perturbation Theory for Linear Operators, (Springer, Berlin) (1966). 

[27] W. D. Heiss, Repulsion of resonance states and exceptional points, Phys. Rev. E 61, 929-932 

( 2000 ). 

[28] 1. Rotter, Exceptional points and double poles of the S matrix, Phys. Rev. E 67, 026204 
(2003). 

[29] U. Giinther, 1. Rotter, and B. F. Samsonov, Projective Hilbert space structures at 
exceptional points. Journal of Physics A: Mathematical and Theoretical 40, 8815 (2007). 

[30] W. D. Heiss and H. Harney, The chirality of exceptional points, The European Physical 
Journal D - Atomic, Molecular, Optical and Plasma Physics 17, 149-151 (2001). 

[31] H. Mehri-Dehnavi and A. Mostafazadeh, Geometric phase for non-Hermitian Hamiltonians 
and its holonomy interpretation, Journal of Mathematical Physics 49, 082105 (2008). 

[32] 1. Rotter, A non-Hermitian Hamilton operator and the physics of open quantum systems, 
Journal of Physics A: Mathematical and Theoretical 42, 153001 (2009). 

[33] W. D. Heiss, The physics of exceptional points. Journal of Physics A: Mathematical and 
Theoretical 45, 444016 (2012). 


- 13 - 




A unifying E2-quasi-exactly solvable model 


[34] W. D. Heiss and G. Wunner, Fano-Feshbach resonances in two-channel scattering around 
exceptional points, The European Physical Journal D 68 (2014) 284. 

[35] J. Favard, Sur les polynomes de Tchebicheff., C. R. Acad. Sci., Paris 200, 2052-2053 (1935). 

[36] F. Finkel, A. Gonzalez-Lopez, and M. A. Rodriguez, Quasiexactly solvable potentials on the 
line and orthogonal polynomials, J. Math. Phys. 37, 3954-3972 (1996). 

[37] A. Krajewska, A. Ushveridze, and Z. Walczak, Bender-Dunne Orthogonal Polynomials 
General Theory, Mod. Phys. Lett. A 12, 1131-1144 (1997). 


- 14 - 




